Scaling Tests of the Cross Section for Deeply Virtual Compton Scattering 
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We present the first measurements of the ep — > epy cross section in the deeply virtual Compton 
scattering (DVCS) regime and the valence quark region. The Q 2 dependence (from 1.5 to 2.3 
GeV 2 ) of the helicity-dependent cross section indicates the twist-2 dominance of DVCS, proving 
that generalized parton distributions (GPDs) are accessible to experiment at moderate Q 2 . The 
helicity-independent cross section is also measured at Q 2 = 2.3 GeV 2 . We present the first model- 
independent measurement of linear combinations of GPDs and GPD integrals up to the twist-3 
approximation. 

PACS numbers: 13.60.Fz, 13.40.Gp, 14.20.Dh, 13.60.Hb 



Measurements of electro-weak form factors determine 
nucleon spatial structure, and deep inelastic scattering 
(DIS) of leptons off the nucleon measures parton dis- 
tribution functions, which determine longitudinal mo- 
mentum distributions. The demonstration by Ji [H, 
Radyushkin and Mueller et al. of a formalism 
to relate the spatial and momentum distributions of the 
partons allows the exciting possibility of determining spa- 
tial distributions of quarks and gluons in the nucleon as 
a function of the parton wavelength. These new struc- 



ture functions, now called generalized parton distribu- 
tions (GPDV, became of experimental interest when it 
was shown [1[ that they are accessible through deeply 
virtual Compton scattering (DVCS) and its interference 
with the Bethe-Heitler (BH) process (Fig. Q]). Figure Q] 
presents our kinematic nomenclature. DVCS is defined 
kinematically by the limit — i -C Q 2 and Q 2 much larger 
than the quark confinement scale. 

The factorization proofs 0, 0] confirmed the connec- 
tion between DVCS and DIS. Diehl et al. @ showed 
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TABLE I: Experimental ep —> ep'y kinematics, for incident beam energy E — 5.75 GeV. 9 q is the central value of the q-vector 
direction. The PbF2 calorimeter was centered on 8 q for each setting. The photon energy for q' || q is E 1 . 
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FIG. 1: Lowest-order QED diagrams for the process ep — > 
ep7, including the DVCS and Bethe-Heitler (BH) amplitudes. 
The external momentum four-vectors are defined on the di- 
agram. The virtual photon momenta are q = k — k' in the 
DVCS- and A = q — q' in the BH-amplitudes. The invariants 
are: W 2 = (q + p) 2 , Q 2 = -q 2 > 0, t = A 2 , x Bi = Q 2 /(2p-q), 
and the DVCS scaling variable £ = —q 2 /(q-P) w a;Bj/(2— TBj), 
with q = (q + q')/2 and P = p + p'. 



that the twist-2 and twist-3 contributions in the DVCS- 
BH interference terms (the first two leading orders 
in l/Q) could be extracted independently from the 
azimuthal-dependence of the helicity-dependent cross 
section. Burkardt 0] showed that the t-dependence of the 
GPDs is the Fourier conjugate to the transverse spatial 
distribution of quarks in the infinite momentum frame as 
a function of momentum fraction. Ralston and Pire 
Diehl [§] and Belitsky et al. [1(3] extended this interpreta- 
tion to the general case of skewness £ ^ 0. The light-cone 
wave function representation by Brodsky et al allows 
GPDs to be interpreted as interference terms of wave 
functions for different parton configurations in a hadron. 

These concepts stimulated an intense experimental ef- 
fort in DVCS. The HI 0,13 and ZEUS [3] collabora- 
tions measured the cross section for xb) ~ 10 -3 . The 
HERMES collaboration measured relative beam helic- 
ity [HI and beam-charge asymmetries 16,E3. Relative 
beam-helicity [l8| and longitudinal target [19( asymme- 



tries were measured at the Thomas Jefferson National 
Accelerator Facility (JLab) by the CLAS collaboration. 

Extracting GPDs from DVCS requires the fundamen- 
tal demonstration that DVCS is well described by the 
twist-2 diagram of Fig. Q] at finite Q 2 . This letter reports 
the first strong evidence of this cornerstone hypothesis, 
necessary to validate all previous and future GPD mea- 
surements using DVCS. We present the determination of 
the cross section of the ep — ► epj reaction for positive and 
negative electron helicity in the kinematics of Table |TJ 



The E00-110 [20| experiment ran in Hall A [2jJ at 
JLab. The 5.75 GeV electron beam was incident on 
a 15 cm liquid H2 target. Our typical luminosity was 
10 37 /cm 2 /s with 76% beam polarization. We detected 
scattered electrons in one high resolution spectrometer 
(HRS). Photons above a 1 GeV energy threshold (and 77 
coincidences from tt° decay) were detected in a 11 x 12 ar- 
ray of 3 x 3 x 18.6 cm 3 PbF2 crystals, whose front face was 
located 110 cm from the target center. We calibrated the 
PbF2 array by coincident elastic H(e, e' Calo pHRs) data. 
With (elastic) k' = 4.2 GeV/c, we obtain a PbF2 resolu- 
tion of 2.4% in energy and 2 mm in transverse position 
(one-cr). The calibration was monitored by reconstruc- 
tion of the 7T° — > 77 mass from H(e, e'TT°)X events. 

We present in Fig. [2] the missing mass squared ob- 
tained for H(e, e'j)X events, with coincident electron- 
photon detection. After subtraction of an accidental co- 
incidence sample, we have the following competing chan- 
nels in addition to H(e, e'j)p: ep — > eir°p, ep — > eir^Nit, 
ep — » ejNir, ep — > ejNmT . . .. From symmetric (lab- 
frame) 7r°-decay, we obtain a high statistics sample of 
H(e, e'ir°)X' events, with two photon clusters in the 
PbF2 calorimeter. From these events, we determine the 
statistical sample of [asymmetric] H(e, e'j)jX' events 
that must be present in our H(e, e'j)X data. The solid 
M\ spectrum displayed in Fig. [2] was obtained after sub- 
tracting this 7T° yield from the total (stars) distribution. 
This is a 14% average subtraction in the exclusive window 
defined by M\ cut in Fig. O Depending on the bin in 
</> 77 and t, this subtraction varies from 6% to 29%. After 
our 7T° subtraction, the only remaining channels, of type 
H(e, e'7)./V7r, Nifn, etc. are kinematically constrained 
to M\ > (M + m w ) 2 . This is the value (Af| cut in 
Fig. [2]) we chose for truncating our integration. Resolu- 
tion effects can cause the inclusive channels to contribute 
below this cut. To evaluate this possible contamination, 
we used an additional proton array (PA) of 100 plastic 
scintillators. The PA subtended a solid angle (relative to 
the nominal direction of the q vector) of 18° < 9 1P < 38° 



and 45° < 



tSyp 



= 180° 



< 315°, arranged in 5 



rings of 20 detectors. For H(e, e'^j)X events near the 
exclusive region, we can predict which block in the PA 
should have a signal from a proton from an exclusive 
H(e, e'jp) event. Open crosses show the X = (p + y) 
missing mass squared distribution for H(e, e'^p)y events 
in the predicted PA block, with a signal above an ef- 
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FIG. 2: (color online). Missing mass squared for H(e, e ,r y)X 
events (stars) at Q 2 = 2.3 GeV 2 and -t £ [0.12,0.4] GeV 2 , 
integrated over the azimuthal angle of the photon <j!> 77 . The 
solid histogram shows the data once the H(e, e'rfjX' events 
have been subtracted. The other histograms are described in 
the text. 



fective threshold 30 MeV. Squares show our inclusive 
yield, obtained by subtracting the normalized triple co- 
incidence yield from the H(e, e /r y)X yield. The dotted 
curve shows our simulated H(e, e'^/)p spectrum, includ- 
ing radiative and resolution effects, normalized to fit the 
data for M\ < M 2 . Triangles show the estimated inclu- 
sive yield obtained by subtracting the simulation from 
the data. Squares and triangles are in good agreement, 
and show that our exclusive yield has less than 3% con- 
tamination from inclusive processes. 

To order twist-3 the DVCS helicity-dependent (<i£) 
and helicity-independent (da) cross sections are [22[ : 



d 4 a 
d 4 * 



d 4 Z(\DVCS\ 2 ) 



sin(0 77 )r I 1 m Im [C 1 ^)] 
sin(2</> 77 )r!, m lm [C 1 '(^ eff )] 



d 4 $ 



(1) 



d 4 $ 

dM\BH\ 2 ] 
d 4 $ 



rf 4 $ 



d 4 a(\DVCS\ 2 ) 



d 4 $ 

„ A Ree [C 1 + AC 1 ] (J 7 ) 



pRe 



+ r^°Ree [C'(F)] - cos(0 77 )rf c Re [C 1 (T)] 
+ cos(20 77 )r^ o Re [C 1 ^]] , (2) 

where d 4 <& = dQ 2 dxBjdtd(f>~ n and the azimutal an- 
gle 77 of the detected photon follows the "Trento- 
Convention" [23| • The rJJ e ' Im are kinematic factors with 
a 77 dependence that arises from the electron propa- 
gators of the BH amplitude. The C 1 and AC 1 angular 
harmonics depend on the interference of the BH ampli- 
tude with the set T = {H, £, Ti., £} of twist- 2 Compton 
form factors (CFFs) or the related set J-° s of effective 



C I (f ee )=F 1 H eS + $G M H 
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[C 1 + AC 1 ] (T) =F 1 H 



AM 2 



F 2 £ - ?Gm [H + £]. (5) 



Fi(t), F 2 (t) and G M (t) = Fi(t) + F 2 (t) are the elastic 
form factors. CFFs are defined in terms of the GPDs 
Hf, Ef, Hf, and Ef, defined for each quark flavor /. 
For example (/ 6 {u, d, s}): 



V I d.r 



2.T 



t)-H f (-t,t,t)] 
H f (x,S,t)\. (6) 



Thus, the DVCS helicity-dependent and helicity- 
independent cross sections provide very distinct and com- 
plementary information on GPDs. On one hand, dS mea- 
sures the imaginary part of the BH-DVCS interference 
terms and provides direct access to GPDs at x = £. On 
the other hand, da determines the real part of the BH- 
DVCS interference terms and measures the integral of 
GPDs over their full domain in x. This real part of the 
BH-DVCS interference term is the same interference term 
that can be obtained by measurements of the difference 
of electron and positron (or fi ± ) DVCS cross sections. 

The twist- 2 and twist-3 CFFs are matrix elements of 
quark-gluon operators and are independent of Q 2 (up to 
logarithmic QCD evolution). Their Q 2 variation mea- 
sures the potential contamination from higher twists. 

The helicity-independent cross section also has a 
cos(30 77 ) twist-2 gluon transversity term. We expect 
this term to be small, and do not include it in our analy- 
sis. We neglect the DVCS 2 terms in our analysis. There- 
fore, our results for Im[C x ] and Re[C x ] may contain, re- 
spectively, twist-3 and twist-2 DVCS 2 terms, which en- 
ter with similar </> 77 dependence. However, the DVCS 2 
terms in both da and dT, are kinematically suppressed 
by at least an order of magnitude in our kinematics [13] , 
because they are not enhanced by the BH amplitude. In 
any case, the terms we neglect do not affect the cross 
sections we extract, which are accurately parametrized, 
within statistics, by the contributions we included. 

Our simulation includes internal bremsstrahlung in the 
scattering process and external bremsstrahlung and ion- 
ization straggling in the target and scattering chamber 
windows. We include spectrometer resolution and ac- 
ceptance effects and a full GEANT3 simulation of the 
detector response to the DVCS photons and protons. 
The spectrometer acceptance is defined for both the data 
and simulation by a i?-function cut p3 |. Radiative cor- 
rections for virtual photons and unresolved real pho- 
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TABLE II: Angular Harmonics fit results, Im and Re parts, and their statistical uncertainties. 

Q 2 \{t) (GeV 2 ) t = -0.33 -0.28 -0.23 -0.17 

Im 1.5 2.1 ±0.3 2.1 ±0.3 2.0 ± 0.2 3.2 ± 0.2 

[C 1 ^)] 1.9 1.9 ±0.2 2.3 ±0.2 2.5 ± 0.2 3.2 ±0.2 

2.3 2.1 ±0.2 2.4 ±0.2 2.6 ± 0.2 3.3 ± 0.3 

Im 1.5 2.8 ±2.0 2.5 ± 2.0 0.1 ± 2.1 0.6 ± 2.4 

[C I (^ rcff )] 1.9 0.3 ±1.4 3.8 ±1.5 -0.9 ±1.8 4.7 ± 2.7 

2.3 5.3 ±1.6 0.7 ±1.8 0.2 ± 2.5 4.0 ±4.6 

Q 2 — 2.3 GeV 2 , Re part of Angular Harmonics 

C(T) -2.4 ±0.1 -2.0 ±0.1 -1.7 ±0.1 -0.7 ±0.2 

[C + AC](F) 0.1 ±0.1 0.8 ±0.1 1.6 ±0.1 2.5 ±0.1 

[COF 80 )] -1.4 ±0.5 0.6 ±0.6 1.0 ±0.8 3.4 ±1.4 



tons are applied according to the VCS (BH±Born am- 
plitude) specific prescriptions of Ref. [25| • This results in 
a global correction factor (independent of 4> 11 or helicity) 
of 0.91 ±0.02 applied to our experimental yields. Within 
the quoted uncertainty, this correction is independent of 
the kinematic setting. 

For each (Q 2 : XB i: t) bin, we fit the Re and 
Im parts (as appropriate) of the harmonics C n E 
{C x (J r ),C :r (J :eS ), [C x + AC 1 ] (J 7 )} as independent pa- 
rameters. In Kin-1 and Kin-2, due to the lower photon 
energy EL (Tabled, our acceptance, trigger, and readout 
did not record a comprehensive set of ep — > eir°X events. 
For those events we were able to reconstruct, we found 
only a few percent contribution to <ffi, but a larger con- 
tribution to da. For Kin-1, 2, we only present results on 
<iE. Our systematic errors in the cross-section measure- 
ments are dominated by the following contributions: 3% 
from HRSxPbF2 acceptance and luminosity; 3% from 
H(e,e'j)jX (tt°) background; 2% from radiative correc- 
tions; and 3% from inclusive H(e, e'^f)NTr . . . background. 
The total, added in quadrature, is 5.6%. The dS results 
contain an additional 2% systematic uncertainty from the 
beam polarization. In order to compute the BH contri- 
bution in the da analysis we used Kelly's parametrization 
of form factors [26| , which reproduce elastic cross-section 
world data in our t range with 1% error and 90% CL. 

For one (Q 2 ,XB},t) bin, Fig. [3] shows the helicity- 
dependent and helicity-independent cross sections, re- 
spectively. We notice that the twist-3 terms make only 
a very small contribution to the cross sections. Note 
also that da is much larger than the BH contribution 
alone, especially from 90° to 270°. This indicates that 
the relative Beam Spin Asymmetry BSA = d A Y>/d A a 
cannot be simply equated to the imaginary part of the 
BH-DVCS interference divided by the BH cross section. 
Table |H] lists the extracted angular harmonics. Fig- 
ure H] (Left) shows the Q 2 dependence of the imaginary 
angular harmonic Imp 1 ] over our full t domain, with 
(t) = -0.25 GeV 2 ((t) varying by ±0.01 GeV 2 over Kin 
1-3). 



The absence of Q 2 dependence of ImfC 1 ^)] within its 
3% statistical uncertainty provides crucial support for the 
dominance of twist-2 in the DVCS amplitude. Indeed, it 
sets an upper limit < 10% to twist-4 and higher contri- 
butions. Im[C x (J-")] is thereby a direct measurement of 
a linear combination of GPDs. The two twist-2 angular 
harmonics extracted from da determine distinct combi- 
nations of GPD integrals, providing most valuable com- 
plementary information on GPDs. As noted above, the 
angular harmonic terms in Table |TT] may include contri- 
butions from kincmatically suppressed bilinear DVCS 2 
terms omitted in our analysis. In our experiment the 
acceptance-averaged ratios of the kinematic coefficients 
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FIG. 3: (color online). Data and fit to d 4 Y 1 /[dQ 2 dxBjdtd(f>~ / ~ / ], 
and d 4 a/[dQ 2 dxBjdtd<f> 11 ], as a function of 4>n- Both are in 
the bin (Q 2 ,t) = (2.3,-0.28) GeV 2 at (x B j> = 0.36. Error 
bars show statistical uncertainties. Solid lines show total fits 
with one-o" statistical error bands. Systematic uncertainty 
is given in the text. The dot-dot-dashed line is the |BH| 2 
contribution to d 4 a. The short-dashed lines in d 4 E and d 4 a 
are the fitted Im and Re parts of C T (T), respectively. The 
long-dashed line is the fitted Re[C z ± AC X ]{T) term. The 
dot-dashed curves are the fitted Im and Re parts of C I (^ reff ). 
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Im C(F) , Im C (F ) (integrated over t) ■ 
lm C'(F) Q 2 =1.5,1.9,2.3 GeV 2 
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FIG. 4: (color online) Left: Q 2 dependence of Im parts of 
(twist-2) C X {T) and (twist-3) C T (^F eS ) angular harmonics, av- 
eraged over t. The horizontal line is the fitted average of 
Im[C (J 7 )]. Right : Extracted real and imaginary parts of 
the twist-2 angular harmonics as functions of t. The VGG 
model curves are described in the text. Note the sign of 
-[C z + AC X ](.F) (data and VGG). Superposed points in both 
panels are offset for visual clarity. Their error bars show sta- 
tistical uncertainties. 



of the DVCS 2 terms to the BH-DVCS terms are below 
1.2% for dT, and below 4.5% for da. The cross-section 
measurements we present are accurate, to the quoted 
uncertainty, and not sensitive within statistics to the ne- 
glected terms in their harmonic analysis. 

Figure 01 (Right) displays the twist-2 angular harmon- 
ics of Table [TT| (Re and Im parts) as functions of t, to- 
gether with the predictions from a model by Vander- 
haeghen, Guichon and Guidal (VGG) [13, H, SH- The 
VGG model (twist-2 contribution only, profile parameter 
b V ai = b sea = 1, Regge parameter a' = 0.8 GeV -2 , GPD 
Ej = 0) is in qualitative agreement with the Im[C z (J? r )] 
data, but significantly under-predicts the principal-value 
integrals (Re parts of the angular harmonics). 

In summary, we present the first explicit demonstra- 
tion of exclusivity in DVCS kinematics. We also present 
the first measurements of DVCS cross section in the va- 
lence quark region. From the Q 2 dependence of the 
angular harmonics of the helicity-dependent cross sec- 
tion, we provide solid evidence of twist-2 dominance in 
DVCS, which makes GPDs accessible to experiment even 
at modest Q 2 . This result supports the striking predic- 
tion of perturbative QCD scaling in DVCS [HQ. As a 
consequence of this evidence for scaling in the exclusive 
channel, and our separate determination of the helicity- 
dependent and helicity-independent cross sections, we ex- 
tract for the first time a model-independent combination 
of GPDs and GPDs integrals. 

We acknowledge essential work of the JLab accelerator 
staff and the Hall A technical staff. This work was sup- 
ported by DOE contract DOE-AC05-06OR23177 under 
which the Jefferson Science Associates, LLC, operates 
the Thomas Jefferson National Accelerator Facility. We 



[1 
[2 
[3 

[4; 

[5 
[6 

ir 

[8 
[9 
[10 
[11 
[12 
[13 
[14 

[is; 

[16 
[17 

[is: 

[19 

[20 
[21 
[22 

[23 

[24 

[25 

[26 
[27 

[28 

[29 



X.-D. Ji, Phys. Rev. Lett. 78, 610 (1997), hep- 
ph/9603249. 

A. V. Radyushkin, Phys. Rev. D56, 5524 (1997), hep- 
ph/9704207. 

D. Mueller, D. Robaschik, B. Geyer, F. M. Dittes, 
and J. Horejsi, Fortschr. Phys. 42, 101 (1994), hep- 
ph/9812448. 

J. C. Collins and A. Freund, Phys. Rev. D59, 074009 
(1999), hep-ph/9801262. 

X.-D. Ji and J. Osborne, Phys. Rev. D58, 094018 (1998), 
hep-ph/9801260. 

M. Diehl, T. Gousset, B. Pire, and J. P. Ralston, Phys. 

Lett. B411, 193 (1997), hep-ph/9706344. 

M. Burkardt, Phys. Rev. D62, 071503 (2000), hep- 

ph/0005108. 

J. P. Ralston and B. Pire, Phys. Rev. D66, 111501 
(2002), hep-ph/0110075. 

M. Diehl, Eur. Phys. J. C25, 223 (2002), hep- 
ph/0205208. 

A. V. Belitsky, X.-D. Ji, and F. Yuan, Phys. Rev. D69, 
074014 (2004), hep-ph/0307383. 

S. J. Brodsky, M. Diehl, and D. S. Hwang, Nucl. Phys. 
B596, 99 (2001), hep-ph/0009254. 

A. Aktas et al. (HI), Eur. Phys. J. C44, 1 (2005), hep- 
ex/0505061. 

C. Adloff et al. (HI), Phys. Lett. B 517, 47 (2001), hep- 
ex/0107005. 

S. Chekanov et al. (ZEUS), Phys. Lett. B573, 46 (2003), 
hep-ex/0305028. 

A. Airapetian et al. (HERMES), Phys. Rev. Lett. 87, 
182001 (2001), hep-ex/0106068. 

F. Ellinghaus (HERMES), Nucl. Phys. A711, 171 (2002), 
hep-ex/0207029. 

A. Airapetian et al. (HERMES) (2006), hep-ex/0605108. 
S. Stepanyan et al. (CLAS), Phys. Rev. Lett. 87, 182002 
(2001), hep-ex/0107043. 

S. Chen et al. (CLAS), Phys. Rev. Lett. 97, 072002 

(200 6), hep-ex/0605012. 

URL http : //hallaweb . j lab. org/experiment /DVCS 
J. Alcorn et al., Nucl. Instrum. Meth. A522, 294 (2004). 
A. V. Belitsky, D. Mueller, and A. Kirchner, Nucl. Phys. 
B629, 323 (2002), hep-ph/0112108. 

A. Bacchetta, U. D'Alesio, M. Diehl, and C. A. Miller, 
Phys. Rev. D70, 117504 (2004), hep-ph/0410050. 
M. Rvachev, Hall A Technical Note JLab-TN-01-055, Jef- 
ferson Lab (2001). 

M. Vanderhaeghen et al., Phys. Rev. C62, 025501 (2000), 
hep-ph/0001100. 

J. J. Kelly, Phys. Rev. C70, 068202 (2004). 

M. Vanderhaeghen, P. A. M. Guichon, and M. Guidal, 

Phys. Rev. D60, 094017 (1999), hep-ph/9905372. 

K. Goeke, M. V. Polyakov, and M. Vanderhaeghen, Prog. 

Part. Nucl. Phys. 47, 401 (2001), hep-ph/0106012. 

M. Guidal, M. V. Polyakov, A. V. Radyushkin, and 

M. Vanderhaeghen, Phys. Rev. D72, 054013 (2005), hep- 

ph/0410251. 



